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Abstract 

We obtain a new lower bound on the radius of the largest ball of separable unnormalized 
states around the identity matrix for a 3-qubit system. This also enables us to improve the 
corresponding lower bounds for multi-qubit systems. These bounds are approximately 5% better 
than the previously known ones. As a by-product, we compute the radius of the largest ball that 
fits into the triple projective tensor product of the unit ball in R 3 . 

1 Introduction 

This work deals with lower bounds on the largest ball of separable unnormalized states around the 
identity matrix for a multi-partite quantum system consisting of qubits. For a detailed introduction 
and a motivation of the paper we refer to [3] . 

In the latest version of [2] L. Gurvits and H. Barnum derived the bounds obtained in [3] by different 
means than those exploited in [3] . Their ideas are applicable also in the case of systems which are more 
complex than multi-qubit systems. In this contribution we improve these bounds, but again our result 
is limited to multi-qubit systems. The result is based on the one hand on the ideas in [3], which apply 
in general to the multi-qubit case, and on the other hand on an idea that is specific to the 3-qubit case. 
The latter will be the main subject of the present paper. As a by-product, we describe those points on 
the boundary of the projective tensor product (see pQ) of 3 unit balls in R 3 which are closest to the 
origin ( Theorems 13 . 813 . 9|l . 

We prove that for a 3-qubit system a ball of radius ^/16/19 around the identity matrix consists only 
of separable elements ^Theorem I3.17|) . as opposed to the best bound -\/4/5 known previously EMU]- 
For systems consisting of more than 3 qubits we obtain an improvement of more than 5% with respect 
to the best bounds known before |2j. Namely, we prove that for an m-qubit system, a ball of radius 
^jx-^ ; 3+1 around the identity matrix consists only of separable elements ( Theorem 14. 

The paper is organized as follows. In the next section we precisely define second order cones, 
positive and separable cones and tensor products of more than two cones. Further we provide some 
simple properties of these objects and state some results from [3], which will be needed later on. In 
Section 3 we consider the cone of separable states for a system of 3 qubits and the related projective 
tensor product of unit balls in R 3 . In the last section we apply the results obtained for 3 qubits to 
study the cone of separable unnormalized states of a multi-qubit system. 



2 Definitions and preliminaries 

In this section we define the Lorentz cone, the cones of positive maps and the separable cones and 
provide some simple properties of these objects. Further, we show how these cones are related to the 
cone of unnormalized separable states for a multi-qubit system. 

2.1 Separable elements and positive maps 

Denote the standard basis elements of R n by eft , . . . , e"_ 1 . Tensor products of real vector spaces will 
have multi- indexed coordinates and basis elements. In particular, e/y, k = 0, . . . , n — 1; I = 0, . . . , m — 1 
will be the canonical basis elements of the space 1" (g) R m . We identify this space with the space of real 
n x m matrices. An element x € R" (g> R m will have coordinates x/-i , k = 0, . . . , n — 1; I = 0, . . . , m — 1, 
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corresponding to the entries of the representing matrix. In this representation a product element 
x ® y £ R™ <E> R m is given by the rank 1 matrix rry T . Tensor products of three real vector spaces will 
have basis elements ejki with three indices, and a point x from such a tensor product space will have 
coordinates Xju- The coordinates of a product x — a ® b ® c will then be given by Xjki — Oj-bftCj. 

Definition 2.1. Let K,K' be regular convex cones, residing in finite-dimensional real vector spaces 
E, E' . We call a linear map M : E -> £" K-to-K' positive if M[Jf] C if'. The set of if-to-if' positive 
maps forms a regular convex cone, the K-to-K' positive cone. 

Definition 2.2. Let K\, . . . ,K n be regular convex cones, residing in finite-dimensional real vector 
spaces Ei, ... , E n . Then an element w £ Ei®- ■ -®E n of the tensor product space is called Ki®- ■ -®K n - 
separable if it can be represented as a finite sum ^fe=i v i ® ' ' ' ® v « 01 product elements such that 
Vj S if; for all fc = 1, . . . , N; I — 1, . . . , n. The K\ ® ■ ■ ■ ® X„-separable elements form a regular convex 
cone, the K\ ® ■ ■ ■ ® K n -separable cone. We denote this cone by K\ ® ■ ■ ■ ® K n . 

It is easily seen that K i ® ■ ■ ■ ® K n = K\ ® (K 2 ■ ■ ■ ® K n ) = {K\ ® ■ ■ ■ ® K n -i) ® 2f„. 

We represent linear maps from R m to R" bynxra matrices. Any linear map M : R m — ► R" can be 
viewed as linear functional on the space R"<g>R' Ti . On the product elements x®y £ R"<g)R m it is defined 
as (M,x ® y) — (M(y),x), where the scalar product on the right hand side is the Euclidean scalar 
product in R™. The vector spaces underlying the cone of positive maps and the separable cone are 
hence dual to each other. With the above matrix representation of the elements of 1" ® R m the scalar 
product of a linear map M : W n -> 1" with an element S £ R™ ® R m is given by (M, S) = tr (M T S) = 
tr (S T M). Here T denotes the transpose. In general we identify the tensor product space R" 1 ®- • -®R ram 
with its dual by means of the Euclidean scalar product (x, y) = 53^=0 ' ' 'Sfe"*=o x ki...k m Vki...k m - 
Remark 2.3. Let E be a real vector space equipped with a scalar product (•,-) and let K C E be a 
convex cone. Then the dual cone K* is defined as the set of elements y £ E such that (x, y) > for 
all x £ K. 

We have the following standard result in the theory of positive maps. 

Theorem 2.4. Let K,K' be regular convex cones. Then the K ® K' -separable cone is dual to the cone 
of K' -to-K* positive maps with respect to the scalar product defined above. 

Proof. We have to show that the K'-to-K* positive maps are precisely those maps which are non- 
negative on all K ® inseparable elements. By definition of the separable cone it suffices to consider 
separable product elements v (g> v' with v £ K, v' £ K' . 

Let M be a positive map and let v £ K, v' £ K' . Then by definition (M, v®v') = (M(v r ), v). But 
v' £ K', and M is K'-to-K* positive. Hence M(v') £ K*. Since v £ K, we then have (M(v'),v) > by 
definition of the dual cone K* . Therefore any positive map is in the dual cone to the K ® K '-separable 
cone. 

Let now M be in the dual cone to the K ® /("'-separable cone. Then (M(V), v) > for all v £ K, 
v' £ K' . But this is equivalent to M(v') £ K* for all v' £ K', hence M is a K'-to-K* positive map. □ 

The dual cone (K\ (g> Ki ® Ks)* to the K\ ® K2 ® ^-separable cone is hence the cone of K^-to-P 
positive maps, where P is the cone of K%-to-K \ positive maps. 

Let A : 1™ — > 1" be a non-degenerate linear map, represented by an n x n matrix. We have the 
following result. 

Lemma 2.5. Let K C R" be a regular convex cone. Then the dual cone to the image A[K] is the 
image A~ T [K*]. 

Proof. Let y £ R™ be an element of the dual space to R n . Then y £ (A[K])* if and only if (y,Ax) = 
y T Ax > for all x £ K. This is in turn equivalent to the condition x T A T y = (x, A T y) > for all 
x £ K, which amounts to the inclusion A T y £ K* . But this inclusion is equivalent to the inclusion 
y£A- T [K*]. □ 

Definition 2.6. Let E\, . . . , E n be real vector spaces and let Aft : E k — > £*, k = 1, . . . , n be linear 
maps. Then the tensor product A\ ® ■ ■ ■ ® A n of these maps is the unique linear map A : E\ ® 
■ ■ ■ ® E n — ► Ei ® ■ ■ ■ ® E n such that for any product element x\ ® ■ ■ ■ ® x n £ E\® ■ ■ ■ ® E n we have 
A(xi ® ■ ■■ ® x n ) = Ai(xi) ® ■ ■■ ® A n (x n ). 

We have the following simple property of tensor products of maps. 

Lemma 2.7. Let K\ C E\, . . . , K n C E n be regular convex cones and let Ak : Ek — > Ek, k = 1, . . . ,n 
be linear maps. Then with A = A\ ® ■ ■ ■ ® A n we have the relation A[K\ ® ■ ■ ■ ® K n ] = A\ [Ki] ® ■ ■ ■ ® 
A n [K n }. □ 

Let us introduce the following notion from p^. 
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Definition 2.8. Let K\ C E\, . . . , K n C E n be convex bodies residing in finite-dimensional real vector 
spaces. Then their projective tensor product K = K\ ® ■ ■ ■ (8 K n C E\® . . . E n is defined as the convex 
hull of the set {x = x\ ® ■ ■ ■ ® x n \ x\ € K\, . . . , x n £ K n }. 

Definition 2.9. Let K S R™ be a set. Then the polar K° of K is defined to be the set {x £ 
l n | |(ar,j/>| < 1} for all y € 

Let <9if denote the boundary of K. The following result follows easily from these definitions. 

Lemma 2.10. Let K be convex, compact and centrally symmetric and define rjf = mm xe QK \x\. Then 
(if ) = K, K° is convex, compact and centrally symmetric and max ie K« \x\ = r^ 1 . Moreover, if y is 
a unit vector, then tk 1 2/ G ?/ and onZj/ if r~^y £ if°. □ 



2.2 Lorentz cones and separable matrices 

Let L n be the n-dimensional standard Lorentz cone, or second order cone, 

L n = |(x ,...,a;„_i) T e R n \x > \j x\ + ■ ■ ■ + x 2 n _^ . 

The Lorentz cone L n is self-dual, i.e. coincides with its dual cone with respect to the standard Euclidean 
metric in R n . By the self-duality of the Lorentz cones a map M : R m — > R™ is L m -to-L n positive if 
and only if M T is L m -to-L n positive. The central ray of L n is generated by the element eg . 

We will denote the L n -to-L n positive cone by P n . Then the dual cone to the triple tensor product 
L n <8> L n (g) L n will be the cone of L n -to-P n positive maps. Therefore an element x is in (L® 3 )* if 
and only if the n x n matrix M(y) defined elementwise by Mij(y) = Xijo + Sfc=i x ijkVk represents 
a L n -to-L n positive map for all vectors y £ R n_1 of norm 1. Clearly for any x £ (X® 3 )* we have 
xooo > 0, and xqqo = implies x = 0. 

Let Ik be the fc x k identity matrix. Denote the space ofnxn hermitian matrices by Tl(n) and the 
cone of positive semidefinite matrices in this space by H + (n). Then we can identify the tensor space 
(H.(n))® m with the space TC(n m ). The tensor product of elements in H(n) amounts to the Kronecker 
product of matrices. The set of unnormalized density matrices for an n-state quantum system is given 
by the cone H + (n). A quantum system consisting of m subsystems with n states each has n m states, 
and the set of its unnormalized density matrices is H + (n m ). The set of unnormalized separable density 
matrices is a subset of H + (n m ), namely the tensor product (7?+(n))® m . 

We will concentrate on the case of multi-qubit systems, i.e. the case n = 2. Note that the space 7i(2) 
can be identified with R via the orthonormal basis of normalized Pauli matrices (and the normalized 
identity matrix I2/V2). Then the cone H + (2) will be identified with the Lorentz cone L4 and the 
cone of separable unnormalized density matrices of an m-qubit system with the cone Lf m . Moreover, 
the identity matrix I2 £ W(2) will be identified with the vector \J~2e\ £ R 4 , and the identity matrix 
Jam = If m £ H(2 m ) with the vector 2" l / 2 e ...o G (R 4 )®" 1 . 

We are interested in the largest radius of a ball around /2" 1 that consists entirely of separable density 
matrices. It is equal to the largest radius of a ball around 2 m / 2 eo...o that is contained in Lf m . 



2.3 Linear images and largest radii of separable balls 

In this subsection we state results from [d] that will be needed further. 

In the sequel we need linear images of Lorentz cones. The Lorentz cone L n can be represented as 
the conic hull of a closed ball in R™ which does not contain the origin. Therefore any linear image of 
L n can be represented as the conic hull of a closed ellipsoid not containing the origin, and any such 
conic hull is the linear image of a Lorentz cone. We refer to such cones as ellipsoidal cones. Moreover, 
by a rotation of R™ such an image can be transformed to some standardized ellipsoidal cone. 

Definition 2.11. Let P be a positive definite real symmetric (n — 1) x (n — l)-matrix. Then the 
standardized ellipsoidal cone generated by P will be 

K eU (P) = [(x ,xi, . . .,x n -i) T I x > Vx T Px, x = (xi, . . . ,x„^i) T | C R". 

It is not hard to establish the following result. 

Lemma 2.12. Let P be a positive definite real symmetric (n— 1) X (n— \)-matrix. Then (K e u(P))* = 
K ell {P- 1 ). ' □ 
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The set of positive definite symmetric (n— 1) x (n— l)-matrices parameterizes the set of standardized 
ellipsoidal cones in M. n . We will need only ellipsoidal cones of a special type, namely conic hulls of balls 
with radius p < 1 around the element These cones are already in the standardized form defined 
above, with the matrix P being a multiple of the identity. The relation between P and p is given by 
the following lemma. 

Lemma 2.13. J3J Lemma 8] The cone generated by a ball of radius p < 1 around the unit vector e$ 
equals the standardized ellipsoidal cone K e u(r~ 2 I n -i) C M" with r = 9 



Let us denote the cone X e/i (r- 2 /„_i) C K n by K£ all (r). Then L n = K^ U {1). Identify R nm with 
M™ l (g)]R rl by identifying the basis vectors e^, k = 0, . . . , m— 1; I = 0, . . . , n— 1, with the orthonormal 
basis of tensor products e™ ® ef. Then the cone K^(r) is generated by a ball centered on e™ <g> e$. 
From [2j we have the following results. 

Lemma 2.14. Corollary 4-] Let n,m £ N, ri,r2 > 0. Let r be the largest number such that the 
inclusion Kffi(r) C K^ u (n) ® tf£ H (r 2 ) folds. Tfon 



= [max{-l + (1 + rl 2 )(1 + r 2 2 ), (min(n,m) - 1)^ 2 rvT 2 }] 



-1/2 



Corollary 2.15. ' L 3, Corollary 6] Let B ri C Tt(m), B r2 C H(n) be balls of radii r\ < \fm, r-i < y/n 
around the corresponding identities L m ,L n and let K\,K2 be the conic hulls of these balls. The largest 
ball around I nm G TL{mn) = 7i(m)<g>7i(n) which is contained in the cone of K\(&Ki-separable matrices 
has radius 



min rif2, — ; . □ 

\ y (min(m 2 , n 2 ) — l)(m — r\)(n — r|) + r\r\ J 



3 Largest ball of separable 3-qubit states 

In this section we obtain a lower bound on the radius of the largest separable ball around the identity 
matrix for a 3-qubit system. This is the main part of the present contribution. 

Lemma 3.1. Let A, B be real symmetric matrices and C a real matrix, all of size nxn. Let Ai, . . . , A n 
be the eigenvalues of A and fii, . . . , p n be the eigenvalues of B, both in decreasing order. Suppose further 
that 

Then \C\% < J2k=i ^k^k- 

Proof. The matrix inequality is fulfilled if and only if C can be represented in the form C = A 1 / 2 MB 1 / 2 , 
where o~ m:ix (M) < 1. Since the 2-norm is a convex function and the extreme points of the set 
{M|<7 max (M) < 1} are the orthogonal matrices, we have \C\% < max^ | uu T =l \A X I 2 XJ B 1 ! 2 ^ and 
|C|| < m'AXjj | uu t-j(A, UBU T ). The set of orthogonal matrices is a smooth manifold, and the sta- 
tionary points of the function (A, U BU T ) on this manifold arc those matrices U that make A and UBU T 
commute. Thus we find an orthogonal matrix V such that A = VAV T and B = VUBU T V T are both 
diagonal. But then the diagonal elements of A are the and the diagonal elements of B are the pk- 
Moreover, we have (A,UBU T ) = (A,B). We obtain \C\\ < maxo. £ s„ J2k=i ^feMo-(fe) = J2k=i ^kPk- 
Here S n is the set of permutations of {1, . . . , n}. □ 

Remark 3.2. The lemma can be extended to the case when A, B are of different size. 

Lemma 3.3. Let M be an n x n matrix with singular values o~\, . . . , a n . Then the matrix 

( aJ n M \ 

{ M T Pin ) 

has spectrum S±S. ± yj + . . . , S±& ± ^2^) 2 + (T 2. 

Proof. Let M — UDV be the singular value decomposition of M. Then we have 

( U T \ / aJ n M \( U \ _ / al n D \ 
{ V ) { M T 0I n ) \ V T ) - \ D [3L n )■ 

Conjugation by an orthogonal matrix does not change the spectrum, and the spectrum of the matrix 
on the right-hand side is exactly 2±i ± y^^^ + o 2 , . . . , 2±£ ± y / ' (^j 2 +0 2 . □ 
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The following lemma is a well-known result in the theory of matrix polynomials, and follows from 
the spectral factorization theorem 0|. 

Lemma 3.4. Let Aq, A±, A2 be real symmetric matrices. Then the condition Ao+cosipAi+simpA 2 ^ 
for all ip G [—71", 7r] is equivalent to the linear matrix inequality 

rX ( A Q + Ai A 2 + X 



Lemma 3.5. Let v = (vi,v 2 ) £K 2 , let Mo, M\, M 2 benxn matrices, such that <7 max (Mo + cos ipM% + 
s'mipM 2 ) < 1 +cos(pv 1 +sin(pv 2 for all <p G [-7r,7r]. Then \M \l + |Mi|jj + \M 2 \ 2 < 2n - (1 - M 2 )(n 
mod 2). 

Proof. Assume the conditions of the lemma. Then we have \v\ < 1. Let ipo be an angle such that v\ — 
\v\ cosipo, v 2 — \v\ sin <p>Q. Then cos(ip + (po)vi + sin(tp + tpo)v 2 = cosip\v\ as trigonometric polynomials in 
ip. Dehne further matrices M[, M 2 by the identity cos((p+(po)Mi+sm((p+ipo)M 2 = cos ipM[ +sin ipM 2 , 
where the left-hand and the right-hand side are also considered as trigonometric polynomials in (p. 
Then we have |Mi|| + |M 2 || = |M{|| + \M 2 \l and a max (M + cos <pM[ +smipM^) < 1 + cos^|w| for all 
tp G [— 7r, 7r]. This is equivalent to the condition 



(1 + cos ip\v\)I n M + cosipM[ + sin<p>M 2 

(Mo + costpM[ + smipM^) T (1 + cosip\v\)I„ 



By Lemma l3.4l this in turn is equivalent to the condition 



3 X x 



-Xl ' , X 2 — —X 2 ,X : 



( (l + MKn 

V mr+x? 



X! 

(l-\v\)In 
(M^) T - X T 

Mf-mr 



hO V p G [-7T,7r]. 



M + M{ 

(1 + MKn 
-X 2 



M2 + X \ 

M - M[ 
X 2 

(l~\v\)In J 



y 0. 



In particular, we have the matrix inequalities 
/ (1 + \v\)I n Xr \ 

\ -x x (i-h)/„;- u ' 

Let Ai , . . . , A n be the singular values of X\ , and fii , 
order. Then we have by Lemmata 13 . II l3~3l 



(i + M)4 

-x 2 



x 2 

(1 - \v\)In 



h 0. (1) 
H n the singular values of X 2 , both in decreasing 



M + M[ M' 2 + X 
X M - M[ 



M' 2 



2(|M 



\M[\ 2 2 



\M' 2 \ 2 2 



mi) 



< 



E 



fe=i 



By ||IJ we have Afc, /ifc < \/l — |v| 2 . If n is even, then we get the bound 



\M \ 2 2 + \M[\l + \M 2 \ 2 2 + \X\ 2 2 <2n. 



If n is odd, then we have A„ = fi n = 0, because a skew-symmetric matrix of odd size is singular. For 
k = 1, . . . , n — 1 we again have the bound Afc, /z/. < yl — |i>| 2 . Hence we obtain 



\M[\ 2 2 



\M0\2 T |J"i| 2 T |i« 2 l2 

The assertion of the lemma now readily follows. 



M'M + \X\ 2 <2n-l 



□ 



Lemma 3.6. Let v G M. m be a column vector and let Mo, Mi, . . . , M m be n x n matrices such that for 
all row vectors x — (x%, . . . , x m ) G K m on the unit sphere we have er max (Mo + X)fc=i x kMk) < 1 + xv. 
Suppose that n is odd. Then \M \\ + Efcli l M fcl2 < 2rl ~ 0- ~ \ v ?) + (m - 2)(n - 1)(1 - |w| 2 ). 
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Proof. Assume the conditions of the lemma. Then we have \v\ < 1. Let U be an orthogonal to x to 
matrix such that Uv = |?;|e™ and define matrices M[ = X^fcLi UikMk and the vector y = (t/i, . . . , y m ) = 
xlI T . Then we have J2T=i l-^fcli = J2T=i l^feli ar± d the condition a maK (M + J2T=i x kMk) < 1 + xv 
is equivalent to CT max (Mo + Sfeli VkM^) < 1 + yi\v\. Let us consider two cases. 

If \v\ = 1, then for y = — (e™) T we have 1 + j/i|i>| = and Mo = M[. Let now fc > 1 and set y = 
(-ef cos 93 + sin ^) T . This yields the condition cr max (M + ^ A^) < f for all tp € [-tt, tt] \ {0}. 

But 1 ^^ ip tends to infinity as tp — > 0, hence M£ = 0. On the other hand, for y = (e™) T we obtain 
cr max (2M ) < 2, and |M || < n. It follows that \M \? 2 + YJk=i \ M 'k\l ^ 2n > which completes the proof 
for this case. 

Let now \v\ < 1. Since n is odd, we have det M = — det(— M) for any n x n matrix. Therefore any 
2-dimensional linear subspace of the space ofnxn matrices has a 1-dimensional subspace consisting of 
singular matrices. We hence find an orthogonal matrix V of size (to — 1) x (m — 1) and n x n matrices 
M 2 , ...,M^ such that Mg +1 = VkiM( +1 for all k = 1, m - 1 and M£, . . . , M^ are singular. 

Define further a vector y' by y[ = y x , (y 2 , ■ • • , y„) = (2/2, • • • , Vm)V T . 

Then we have 0- max (M o + V\M' 1 + J2T=2 Vk^'k) — 1 + 2/1 M f° r au 2/' on the unit sphere. Consider 
the following linear fractional automorphism of the sphere iS™ -1 . 



y'(VT=Wi m + (1 - VT^W)eT(eT) T ) + (eT) T \v\ 
i + y'M 

, = z(y/T=Wl m + ( I - ^T^W)e?(e?) T ) - (e?) T \v\ 



Then the condition o- max (M +y' 1 M[ + J2T=2 Vk^'k) — 1+2/iM 1S equivalent to the condition <r max (5o + 
I2Z=i z kB k ) < 1, where 

B ° = T^W {Mo ~ HM{) ' Bl = — l ~W [M[ HMo) ' Bk = 7i=pF M ^ fc > L 

We have in particular <7 max (£?o ± Bk) < 1 and hence cr max i?fe < f for all fc. Therefore |Bfc|a < ra — 1 
for fc > 2, because at least one singular value of Bk, fc > 2, is zero. This yields the bound < 
(1 - H 2 )(n - 1) for fc > 2. In addition, by the preceding lemma we have \M \ 2 + \M[\\ + \M 2 \\ < 
2n-(l-\v\ 2 ). 

This finally yields |Af |i + EZi \ M k\ 2 2 = \M \ 2 2 + \M[\ 2 + £™ 2 \ M k\l < 2n - (1 - \v\ 2 ) + (to - 
2)(« - 1)(1 - |^| 2 ). □ 

Corollary 3.7. Let v E M 3 be a column vector and let Mo, Mi, M2, M3 be 3 x 3 matrices such that 
Cmax(Mo + X^feli^fcMc) < 1 + xv for all row vectors x on the unit sphere. Then \Mo\ 2 + \Mi\ 2 + 
\M 2 \ 2 2 + \M 3 \ 2 + \v\ 2 <7. " □ 

This bound is exact, as the example 
/ 1 

M = 0, v = 0, Mi = -1 I , M 2 = I -1 I , Ms = I I (2) 






shows. This example together with the preceding corollary yields the following result. 

Theorem 3.8. Let B C R 3 be the unit ball. Then the largest ball centered on the origin in (R 3 )® 3 
that is contained in the projective tensor product _B® 3 = B (g> B B has radius . 

Proof. In this proof we will not deal with cones, but with bounded convex bodies. We hence adapt 
throughout the proof a numeration of the indices that commences with 1. The coordinates of a vector 
x € K 3 will hence be xx, x 2 ,xz- 

Note that B® 3 is convex, compact, and centrally symmetric. A point w S (R 3 )® 3 is in the polar 
of -B® 3 if and only if • k i=1 WjkiXjykZi < 1 for any triple of unit length vectors x,y,z £ R 3 . Define 
3x3 matrices Mj, j — 1,2,3, by (Mj)ki = Wjki- Then the condition w S (B® 3 )° is equivalent to the 
condition cr max (5^ fe=1 XjMj) < 1 for all x G B. We have \w\ 2 — I-^j'IIj hence by the preceding 

corollary and in view of example @ we have max^grg® 3 ) \w\ — \fl . An application of Lemma |2. 101 
concludes the proof. □ 

The following result shows that the matrices Mi, M2, M3 in @ define essentially the only element 
in (£?® 3 )° which has norm \fl . 
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Theorem 3.9. Let B C R 3 be the unit ball and let x e dB® 3 be such that \x\ = \/l/7. Then there 
exist orthogonal 3x3 matrices U, V, W such that x' = 7(U ® V ® W)x is given componentwise by 
x 'jki = (Mi)jk, jjk,l = 1,2,3, where the matrices Mi, M 2 , M 3 are given by Hj)). 

Proof. By Lemma f2. 101 we have only to show that any element y S (_B® 3 )° of norm V7 can be brought 
to the form defined by the matrices Mi, M2, M 3 by a transformation of the form U <g> V (E> VF, where 
£7, V, IF are orthogonal. 

Let y € (£?® 3 )°, |y| 2 = 7 and define matrices Yi,Y 2 ,Y 3 by (il)jfc = yjki, j,k,l = 1,2,3. Assume 
without loss of generality that Y 3 is singular (otherwise we apply an appropriate transformation of the 
form Y 2 1— > co&tpYi + sintpY 3 , Y 3 1— ► — sin + cos <pY 3 , which amounts to the application of a map 
Id ® Id ® W to ?/, where W is a rotation in the (2, 3)-plane). Then \Y 3 \\ < 2 and by Lemma 13.51 
\Yi\l + \Y2\2 < 5. It follows that both inequalities are in fact equalities. In particular, the pendants of 
the matrices M$,X in the proof of Lemma f3 . 51 must be zero, and the pendants of Xi,X 2 must have 
2-norm y/2. This leads to the following condition on Y\, Y2. 



3 Xi 



Xn 



-xi 



{ 


h 


Xi 


Yi 


Yi 


\ 




-x x 


h 


Y 2 


-Yi 




Yi T 




h 


x 2 








~Yi T 


-x 2 


h 


) 



h0; 



It can be written in a more compact form as 

3 Xi = —Xj,X 2 = —X. 



L 2 



h 

r-iY? 



iXi Yi 



~iY 2 
iX 2 



t 0; 



\Xi\ 2 2 = \X 2 \ 2 = 2. 



L 2|2 



There exist orthogonal matrices U, V such that UXiU T and VX 2 V T equal the matrix 




By conjugating above matrix inequality with diag(U, V), we get the matrix inequality 



h + iJ 
V(Yi T - iY 2 T )U T 



U(Yi^ 
h 



iY 2 )V T 
f iJ 



>- 0. 



Then the relation \U{Y X + iY 2 )V T \\ = 5 can hold only if VY X V T = I 3 , UY 2 V T = J. 

Note that cr max (cos tpY\ + sin tpY 3 ) = cr max (cos tpl 3 + sm<pUY 3 V T ) < 1 for all tp. This is possible 
only if UY 3 V T is skew-symmetric. Moreover, we have a- max (cos ipY 2 + sin (pY 3 ) = a max (cos tpUY 2 V T + 
shnpUY 3 V T ) < 1 for all tp. Since both matrices UY 2 V T and UY 3 V T are skew-symmetric, this is 
possible only if they are orthogonal to each other. But then we find another orthogonal matrix U' such 
that 



U'UYiV T {U') T = I 3 , U'UY 2 V T {U') T = 




U'UY 3 V T (U'Y = 




The theorem now follows by multiplying all these three products by the orthogonal matrix diag(\, —1, — 1) 
from the right. □ 



Proposition 3.10. Let x e (Lf )* with xqqq = F Then the inequality 

EL-.fe=i x lk < 7 holds - 



k=l ^00k 



E3 2 
k,l = l x kW 



Proof. Let x satisfy the assumptions of the proposition and define the mapping D : M 4 — > K 4 by 
D = diag(l, —1, —1, —1). Note that the Lorentz cone L4 is invariant with respect to this transformation. 



Then by Lemma |2~7I the cone Lf 3 is invariant with respect to the map D ® D 
it therefore follows that (Lf 3 )* is invariant with respect to the map (D ® D C> 
Define a point x' = (D ® D ® Id)(x). It is given coordinatewise by 



) Id. From Lemma l2~5l 

Id)- T = D(g> D(g> Id. 



l ijk 




i = 3 = 0, 
*¥> 0, j ? 0, 
ij = 0, i + 



and is also in (Lf )*. By convexity of this cone we have y = ^4p- S (Lf ) 
i ^ and all fc; and yojfc = for all j ^ and all k. For all other elements we have yijk 



Note that ?/iofc = for all 

•Eijk • 



7 



Let us now define the vector v € R 3 by Vk = a^oofc, the 3x3 matrix M$ by (Mq)^ = Xyo an d the 
3x3 matrices Mk, k = 1,2,3, by (Mk)ij — Xijk- Then the condition y € (Lf 3 )* is equivalent to the 
condition 

l \ zv Mo + eL**^ ) eP4 Vz = (zi,Z2,Z3)rt3 ' W = L 

This in turn is equivalent to the condition <7 max (Mo + EfcLi z kMk) < 1 + zv for all row vectors z£l 3 , 
I&I = 1. Corollary 13 . 71 then completes the proof. □ 

Consider the inequality Efc=i x ook + Efc,z=i 4o + Eij,fc=i x %k - ^ ^ om tne preceding proposition. 
If we permute the indices and combine all resulting inequalities, then we obtain the following corollary. 

Corollary 3.11. Let x 6 (Lf 3 )* with £ooo = 1- Then the inequality |Efc=i(4ofc + x oko + 4oo) + 

I T,l,i=i(4w + x toi + x Im ) + £l, fc =i 4* - 7 D 

This means that (Lf 3 )* is contained in the ellipsoidal cone defined by this inequality. 

We now make repeatedly use of Lemma 12.141 Consider the space (R 4 )® 2 = R 16 and the cone 
L4, ® £4 = K^ all {l) <g> Ki all (l) in this space. By Lemma \'2. 141 this cone contains the cone K^ u (l/y/3). 
Again by Lemma [2.141 the cone K^ u (l / \/3) ® K^ all (^5/3) contains the cone K^ u (y/5/27). Hence 

Ktiu{\fWx) <^u®u® Kt aU {^m. 

Now note that the cone K^ a n(W5/3) is a linear image of L4. Namely, if D is a linear map given 
by the diagonal matrix diag(l, \ /I, \ /I), then £>[AT 4 qH (-\/5/3)] = L4. Therefore (Id (3 Id ® 



5' V 5' V 5 ^ 

J D)[if b 6 ( 4 H (^5/27)] C (Id® Id® D)[L 4 <g> L 4 <g> K^y/E/S)] = Lf 3 . Let us formalize this result in the 
following lemma. 

Lemma 3.12. Lety 6 (R 4 )® 3 be such thaty^ = 1 and the inequality ^ 3 k=1 (y^ m +yQ k0 )+Y^ 3 k l=1 y k i + 

! (eLi 2/oo fc + eL=i(?/L + vlki) + r:.,.,, , 4- fc ) ^ & ^ T/i ™ » e L f 3 - D 

By Lemma f2. 121 we obtain the following corollary. 

Corollary 3.13. Let x e (Lf 3 )* with x 000 = 1. Then the inequality Efe=i( :c fcoo+ a; ofeo)+EL=i 4o + 

f (ELl *00 fe + EiUlOSfcM + ^) + El 3 ,k=l X lk) < f ^i*. □ 
Again by a permutation of the indices and combining the inequalities we obtain the following 
corollary. 



fcOO ~r ^QfeO ~r~ ^OOkJ 



Corollary 3.14. Let x G (Lf 3 )* with xqqq = 1. XTien i/ie inequality ^Y^k = i( x 

11 eL =1 (4*o + 4 0i + x lki) + 9EL,fc=i 4* ^ 81 /ioMs - D 

Adding 3 times the inequality from Corollary 13 . 1 II we obtain the following theorem. 

Proposition 3.15. Let x G (Lf 3 )* with xooo = !• Then the inequality 7E/!=i(4oo + x oko + -^oofc) + 
6EL=i(4o + x loi + x lki) + 6£<j,*=i 4 fc ^ 5i Wds - D 

Corollary 3.16. TTie inclusions (Lf 3 )* C K^ u (y^V7/2) and K^ u (y / 2/17) C Lf 3 ZioZd. □ 

Note that the conic hull of a ball of radius p around e$ A and the conic hull of a ball of radius Ip 
around le^ 1 , I > 0, are equal. The following theorem is then a consequence of the preceding corollary 
Lemma \2 . 1 31 and the fact that the identity Is € H(8) = R 64 has norm \/~8. 

Theorem 3.17. Let B C H(8) = (7i(2))® 3 be an Euclidean ball of radius ^/j"§ around the identity 
matrix Is- Then B consists of tripartite separable elements. □ 
The theorem yields an improved lower bound on the largest radius of a tripartite separable ball 
around the identity matrix. The best previously known bound was , obtained independently by 
Leonid Gurvits by a different line of reasoning and by the author by using iteratively Corollary 12. 151 
on a triple of Lorentz cones L4. 



4 Largest ball of separable n-qubit states 

Applying Corollary 12.151 further, one gets improved lower bounds on the largest radius of a ball of 
separable multi-qubit states around the identity. 

Following the line of reasoning in Section 5] , we define a sequence pk recursively by 

r>- 2 - — ■ n - 2 --n- 2 -2- fe+1 k > 3 
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We have the explicit expression 



^ a =S(T)* +!r *' pk= i 2fc/2 ' 

54 ./i73fc-3 + i 



2 

By the repeated use of Corollary 12 . 1 51 we obtain the following theorem (see also [3]). 

Theorem 4.1. For k > 3 the expression — -j==== is a lower bound on the radius of the largest 

separable ball of unnormalized multi-partite mixed states of a k-qubit system around the identity matrix 
m the space H(2)® k . □ 

For n-qubit systems with n > 3 we get an improvement of over 5% with respect to the results in 
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